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We introduce a general deterministic model for Apollonian Networks in an iterative fashion. 
The networks have small-world effect and scale-free topology. We calculate the exact results for 
the degree exponent, the clustering coefficient and the diameter. The major points of our results 
indicate that (a) the degree exponent can be adjusted in a wide range, (b) the clustering coefficient 
of each individual vertex is inversely proportional to its degree and the average clustering coefficient 
of all vertices approaches to a nonzero value in the infinite network order, and (c) the diameter 
grows logarithmically with the number of network vertices. 
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I. INTRODUCTION 

In the past few years we have witnessed an upsurge in 
the research of a wide range of complex systems which 
can be described in terms of networks vertices con- 
nected together by edges 0, 0> 01 0] • Various empirical 
studies have shown that most real-life systems exhibit the 
following three properties: power-law degree distribution 
, high clustering coefficient and small average path 
length (APL) which implies that their average interver- 
tex distances grow logarithmically with the number of 
vertices or slower. 

In order to mimic real-world systems, a wide variety of 
models have been proposed [l], 0> L4 > most of which are 
random. Randomness may be in line with the major fea- 
tures of real-life networks, but as mentioned by Barabasi 
et al., it makes harder to gain a visual understanding of 
how networks are shaped, and how do different vertices 
relate to each other pj. Therefore, it would be of ma- 
jor theoretical interest to construct deterministic mod- 
els having similar structural characteristics as systems in 
nature and society. A strong advantage of determinis- 
tic networks is that it is often possible to compute ana- 
lytically their properties, which may be compared with 
empirical data from real-life systems. In the last few 
years, deterministic networks have been intensively stud- 

iedgiinammmnininina. 

Recently, In relation to the problem of Apollonian 
packing, Andrade et al. introduced Apollonian net- 
works [l8j which were also proposed by Doye and Massen 
in |19| . Apollonian networks are part of a deterministic 
growing type of networks and have received much at- 
tention. Except their further properties l2fl. stochastic 
versions plll22| and potential implications |l9ll23| . many 
peculiar results about some famous models upon Apollo- 
nian networks have been found 0, |3 HE @ • 

In this paper, we do an extensive study on Apollo- 
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nian networks. In an iterative way we propose a general 
model for Apollonian networks by including a parame- 
ter. Apollonian networks |18L fl9L f20j are special cases of 
the present model. The deterministic construction of our 
model enables one to obtain the analytic solution for the 
degree distribution, clustering coefficient and the diame- 
ter. Moreover, by tuning the parameter, we can obtain 
scale-free networks with a variety of degree exponents 
and clustering coefficient. 



II. BRIEF INTRODUCTION TO APOLLONIAN 
NETWORKS 

To define Apollonian networks, we first introduce 
Apollonian packing (see Fig. Q]for the case of two dimen- 
sion) . The classic two-dimensional Apollonian packing is 
constructed as follows. Initially three mutually touching 
disks are inscribed inside a circular space which is to be 
filled. The interstices of the initial disks and circle are 
curvilinear triangle to be filled. This initial configuration 
is called generation t — 0. Then in the first generation 
t = 1, four disks are inscribed, each touching all the 
sides of the corresponding curvilinear triangle. For next 
generations we simply continue by setting disks in the 
newly generated curvilinear triangles. The process is re- 
peated indefinitely for all the new curvilinear triangles. 
In the limit of infinite generations we obtain the Apol- 
lonian packing, in which the circular space is completely 
filled with disks of various sizes. The translation from 
Apollonian packing construction to Apollonian network 
generation is quite straightforward: vertices of the net- 
work represent disks and two vertices are connected if 
the corresponding disks are tangent 0, 0] . 

The two-dimensional Apollonian network can be eas- 
ily generalized to high-dimensions (d-dimensional, d > 2) 
2jJl2(j associated with other self-similar packings ]27|. 
The d-dimensional Apollonian packings start with d + 1 
mutually touching (i-dimensional hyperspheres that is en- 
closed within and touching a larger d-dimensional hyper- 
spheres, with d + 2 curvilinear d-dimensional simplex (d- 
simplex) as their interstices, which are to be filled in sue- 
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FIG. 1: (Color online) A two-dimensional Apollonian packing 
of disks. 



cessive generations. If each d-hypersphere corresponds to 
a vertex and vertices are connected if the corresponding 
d-hyperspheres are in contact, then d-dimensional Apol- 
lonian networks are gotten. 

According to the construction process of d-dimensional 
Apollonian packings, in Ref. [20], an generation algo- 
rithm for d-dimensional Apollonian networks was pro- 
posed and further properties of the networks were in- 
vestigated. The d-dimensional Apollonian networks are 
scale- free and display small- world effect 0, |2(]] . The 
degree exponent is 7 = 1 + x^r^ ■ Both their diam- 
eter and average path length grow logarithmically with 
the number of network vertices. The clustering coeffi- 
cient is large, for two- and three-dimensional Apollonian 
networks, it approaches 0.8284 and 0.8852, respectively. 



III. ITERATIVE ALGORITHM FOR THE 
GENERAL MODEL OF APOLLONIAN 
NETWORKS 

Our model is constructed in an iterative fashion. Be- 
fore introducing the algorithm we give the following def- 
initions on a graph (network). The term size refers to 
the number of edges in a graph. The number of vertices 
in the graph is called its order. When two vertices of a 
graph are connected by an edge, these vertices are said 
to be adjacent, and the edge is said to join them. A com- 
plete graph is a graph in which all vertices are adjacent to 
one another. Thus, in a complete graph, every possible 
edge is present. The complete graph with d vertices is 
denoted as Kd (also referred in the literature as d- clique; 
see psjp. Two graphs are isomorphic when the vertices 
of one can be relabeled to match the vertices of the other 
in a way that preserves adjacency. So all d-cliques are 
isomorphic to one another. 

The general model for Apollonian networks after t gen- 
erations are denoted by A(d,t), d > 2,t > 0. Then at 



step t, the considered networks are constructed as follows: 
For t = 0, A(d, 0) is a complete graph Kd+2 (or (d + 2)- 
clique). For t > 1, A(d,t) is obtained from A(d, t — 1). 
For each of the existing subgraphs of A(d,t — 1) that is 
isomorphic to a (d + l)-clique and created at step t — 1, 
m new vertices are created, and each is connected to 
all the vertices of this subgraph. The growing process 
is repeated until the network reaches the desired order. 
When m = 1, the networks reduce to the deterministic 
Apollonian networks [H ITl I20] . 

Now we compute the order and size of the networks. 
Let n v (t), n e {t) and K^ + i) t be the numbers of vertices, 
edges and (d + l)-cliques created at step t, respectively. 
Note that the addition of each new vertex leads to d + 1 
new (d+l)-cliques and d+1 new edges, so we have n v (t) — 
mK(d+\),t-\, n e (t) = {d+ l)n v (t), and K( d +i),t = (d + 
l)n v (t). Thus one can easily obtain -ftT^+ij.t = m(d + 

1) #(d+i),t-i = (d + 2)[m(d + 1)]* (t > 0), n v {t) = (d + 

2) m t (d + I) 1 - 1 (t > 0) and n e (t) = (d + 2)m t (d + if 
(t > 0). So the number of network vertices increases 
exponentially with time, which is similar to many real- 
life networks such as the World Wide Web. From above 
results, we can easily compute the size and order of the 
networks. The total number of vertices N t and edges E t 
present at step t is 

t 

N t = ^n v {U) 

ti=0 

= ^+2)K(d+l)«-l] 
m(d+ 1) - 1 

and 

t 

E t = y^ j n e (U) 

ti=0 

_ m(d + 2)(d+l)[m t (d+l) t - 1] (d + 2)(d+l) 

m(d + l)-l + 2 { ' 

respectively. So for large t, The average degree k t = 
is approximately 2(d + 1). 

IV. TOPOLOGY PROPERTIES OF THE 
NETWORKS 

A. Degree distribution 

When a new vertex i is added to the graph at step tj, 
it has degree d+1 and forms d+1 new (d + l)-cliques. 
Let Kfd+ijii, t) be the number of newly-created (d + 1)- 
cliques at step t with i as one vertex of them, which 
will create new vertices connected to the vertex i at step 
t+1. At step U, K(d + i)(i, U) =d+l. From the iterative 
process, we can see that each new neighbor of i generates 
d new id + l)-cliques with i as one vertex of them. Let 
ki(t) be the degree of i at step t. It is not difficult to find 
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following relations: 

Ah(t) = ki(t) - ki(t - 1) = mK (d+1) (i,t - 1) 

and 

K (d+1) (i,t) = dmK (d+1) (i,t-l) = {d+l)(d m y- u 

Then the degree of vertex i becomes 

t-i 

fo(t) = kjjtj) + m K( d+1) (i, t) 

T=ti 

m(d + l)[(md)*-*' 
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Since the degree of each vertex has been obtained 
explicitly as in Eq. ©, we can get the degree distri- 
bution via its cumulative distribution, i.e. P cum (fc) = 
T,k'>k N ( k '^)/ N t ~ fcl " 7 : where N(k',t) denotes the 
number of vertices with degree k' . The analytic compu- 
tation details are given as follows. For a degree k 



m(d+ l)[(md) 



t- P 



1 



there are n v (p) — (d + 2)m p (d + vertices with this 

exact degree, all of which were born at step p. All vertices 
with birth time at p or earlier have this and a higher 
degree. So we have 



k'>k 



s=0 



m(d + 2)[mP(d+l) p -1] 
m(d+ 1) - 1 



d + 2 



(4) 



As the total number of vertices at step t is given in Eq. 
we have 



m(d + l)[(md) 



t-P 



1 



1-7 



d-1 

m(d+2)[m p (d+l) p -1] , , , ^ 
mid+X)— 1 i - 



t(rf+l)- 

m(d+2)[m t (d+l) t -l] 
m(rf+l)-l 



d + 2 



Therefore, for large t we obtain 

[(md)^] 1 " 7 = Md+l)]^' 

and 

ln[m(d+l)] 



7 ~ 1 + 



ln(md) 



(5) 



For m = 1, Eq. recovers the results previously ob- 
tained in Refs. [TaTll9L l20| . 



FIG. 2: (Color online) The dependence relation of C on d 
and m. 



B. Clustering coefficient 

We can go beyond the degree distribution and obtain 
the analytical expression for clustering coefficient C(k) 
of an individual vertex as a function of its degree k. By 
definition, clustering coefficient Ci of a vertex i is 
the ratio of the total number et of edges that actually 
exist between all ki its nearest neighbors and the number 
ki(ki — l)/2 of all possible edges between them, i.e. Ci = 
2ei/[ki(ki — 1)]. The clustering coefficient of the whole 
network is the average of C^s over all the vertices. 

When a vertex is generated it makes connections to all 
the vertices of a (d + l)-clique whose vertices are com- 
pletely interconnected. It follows that its degree and clus- 
tering coefficient are d + 1 and 1, respectively. In the 
following iterative steps, if its degree increases one by a 
newly created vertex connecting to it, then there must 
be d existing neighbors of it attaching to the new vertex 
at the same time. Thus for a vertex of degree k, we have 



C(k) 



d(d+l) 
2 



+ d(k-d-l) 2d{k 



d+l \ 
2 > 



fc(fc-l) 
2 



k(k - 1) 



(6) 



where the last expression in Eq. JBJ is obtained after 
some algebraic manipulations. From Eq. 10, one can 
easily see that C(k) depends on degree k and d. The 
asymptotic behavior for large k is C(k) ~ 1/fc, which 
implies that C(k) is inversely proportional to the degree. 
Interestingly, we recover in our model the same sca ling 
behavior of C(k) found in other models 0, 0, 
13 IH 13 13 Ifl EllH and some real-life systems [l2|. 

Using Eq. 10, we can obtain the clustering C of the 
networks at step t: 



C 



1 

N't 



t 

E 

r=0 



2d{D r - ^)n v {r) 
D r (D r - 1) 



(7) 
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where the sum is the total of clustering coefficient for all 
vertices and D r = m ( d + 1 )[( m ^~ - 1 shown by Eq. 

|jSJ is the degree of the vertices created at step r. 

In the infinite network order limit (N t — > oo), Eq. (J7J 
converges to a nonzero value. When d = 2, for m = 1, 2 
and 3, C equal to 0.8284, 0.8602 and 0.8972, respectively. 
When m = 2, for d = 2, 3 and 4, C are 0.8602, 0.9017 and 
0.9244, respectively. Therefore, the clustering coefficient 
of our networks is very high. Moreover, similarly to the 
degree exponent 7, clustering coefficient C is determined 
by d and to. Fig. [2 shows the dependence of C on d and 
m. From Fig. 0(a) and (b), one can see that for any fixed 
m, C increases with d. But the dependence relation of C 
on to (see Fig. H](b)) is more complex: (i) when to < 2 
and d < 6, for the same d, C increases with to; (ii) when 
to < 2 and d > 6, for the same d, C decreases with to; 
(iii) when to > 3, for arbitrary fixed d, C increases with 
to. The reason for this complicated dependence relation 
would need further research. 



C. Diameter 

The diameter of a network characterizes the maximum 
communication delay in the network and is defined as 
the longest shortest path between all pairs of vertices. 
In what follows, the notations \x] and |_^J express the 
integers obtained by rounding x to the nearest integers 
towards infinity and minus infinity, respectively. Now we 
compute the diameter of A(d,t), denoted diam(A(d,t)) 
for d > 2 : 

Step 0. The diameter is 1. 

Steps 1 to |~|~|. In this case, the diameter is 2, since 
any new vertex is by construction connected to a (d + 1)- 
clique, and since any (d + l)-clique during those steps 
contains at least i +2 (d even) or ^±1 + 1 odd) vertices 
from the initial (d + 2)-clique A(d, 0) obtained after step 
0. Hence, any two newly added vertices u and v will 
be connected respectively to sets S u and S v , with S u C 
V{A(d,0)) and S v C V{A(d,0)), where V{A(d,0)) is the 
vertex set of (A(d, 0); however, since \S U > |+2 (deven) 
and \S V \ > + 1 (d odd), where \S\ denote thes number 
of elements in set S, we conclude that S u Pi S v ^ 0, and 
thus the diameter is 2. 

Steps [fl + l to d+1. In any of those steps, some newly 
added vertices might not share a neighbor in the original 
(d + 2)-clique A(d, 0) obtained after step 0; however, any 
newly added vertex is connected to at least one vertex of 
the initial clique A(d, 0). Thus, the diameter equals to 3. 

Further steps. Clearly, at each step t > d + 2, the di- 
ameter always lies between a pair of vertices that have 
just been created at this step. We will call the newly 
created vertices "outer" vertices. At any step t > d + 2, 
we note that an outer vertex cannot be connected with 
two or more vertices that were created during the same 
step < t' < t — 1. Moreover, by construction no two 
vertices that were created during a given step are neigh- 
bors, thus they cannot be part of the same (d+ l)-clique. 



Thus, for any step t > d + 2, some outer vertices are con- 
nected with vertices that appeared at pairwise different 
steps. Thus, there exists an outer vertex v t created at 
step t, which is connected to vertices v'^s, 1 < i < t — 1, 
all of which are pairwise distinct. We conclude that vt 
is necessarily connected to a vertex that was created at 
a step to < t — d — I. If we repeat this argument, then 
we obtain an upper bound on the distance from v t to 
the initial clique A(d, 0). Let t = a(d + 1) + /3, where 
1 < (3 < d+1. Then, we see that v t is at distance at most 
a + 1 from a vertex in A(d, 0). Hence any two vertices 
Vt and wt in A(d, t) lie at distance at most 2(a + 1) + 1; 
however, depending on /?, this distance can be reduced 
by 1, since when (3 < |"|] , we know that two vertices cre- 
ated at step j3 share at least a neighbor in A(d, 0). Thus, 
when 1 < p < |"|~|, diam(A(d, t)) < 2(a + 1), while when 
rfl + 1 < P < d+ 1, diam{A{d,t)) < 2(a+ 1) + 1. One 
can see that these distance bounds can be reached by 
pairs of outer vertices created at step t. More precisely, 
those two vertices vt and w t share the property that they 
are connected to d vertices that appeared respectively at 
steps t — l,i — 2, . . .t — d — 1. 

Based on the above arguments, one can easily see that 
for t > d + 2, the diameter increases by 2 every d+1 
steps. More precisely, we have the following result, for 
any d > 1 and t > 1 (when t = 0, the diameter is clearly 
equal to 1): 

diam{A{d,t)) = 2(L^J + 1) + f(d,t) 
d+1 

where f(d, t) = if t - L^jJ (d + 1) < |"f ] , and 1 other- 
wise. 

In the limit of large t, diam(A(d,t)) ~ -£p[, while 
Nt ~ [m(d + 1)]*, thus the diameter is small and scales 
logarithmically with the network order. 

V. CONCLUSION AND DISCUSSION 

In summary, we have proposed and studied a network 
model, which is built in an iterative fashion. At each time 
step, each already existing (d + l)-clique, that is created 
at last time step, produces to new vertices. The iterative 
process leads to a serial of networks with Apollonian net- 
works corresponding to the particular case to = 1. Wc 
have obtained the analytical result for degree exponent, 
clustering coefficient and diameter. The degree exponent 
and the clustering coefficient may be adjusted to various 
values by tuning the parameter to. Moreover, our net- 
works consist of complete graphs, they may represent a 
variety of real-life systems such as movie actor collabora- 
tion networks, scientific collaboration networks and net- 
works of company directors, all of which are composed of 
cliques QHlES- 
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